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Field  Reversals  of  “Paleomagnetic”  Type  in  Coupled  Disk  Dynamos 

John  H.  Mathews  and  W.  K.  Gardner 
Energy  Conversion  Branch 
Electronics  Division 

Rikitake’s  adaptation  of  Bullard's  simplified  disk-dynamo  model  for  the  dynamo  theory  of  the 
geomagnetic  dipole  was  investigated  further  by  means  of  an  analog  computer.  A  modified  model 
having  generalized  interaction,  which  was  suggested  by  an  axial  configuration,  was  also  studied.  A 
striking  similarity  is  observed  betw^n  the  computer  solutions  and  the  geomagnetic  polarity  reversals 
inferred  from  paleomagnetic  studies.  Thus  the  magnetofluiddynamic  origin  of  the  earth’s  field  becomes 
more  plausible. 


INTRODUCTION 

Many  theories  have  been  advanced  to  explain  the 
origin  of  the  earth’s  main  diftole  field.  Some  magneto¬ 
fluiddynamic  process  operating  as  a  homogeneous 
dynamo  in  the  earth’s  core  is  now  generally  accepted, 
with  fair  certainty,  as  the  basic  mechanism  for  a  valid 
theory.  The  formidable  nature  of  the  mathematics  of  the 
problem  has  made  development  of  a  theory  most  diffi¬ 
cult.  As  a  result,  much  of  the  earlier  effort  had  been  on 
steady-state  considerations  of  the  theoretical  homo¬ 
geneous  dynamo.  Some  important  results  are  concerned 
with  the  fundamental  question  of  the  ability  of  the 
steady-state  dynamo  to  maintain  itself.  Backus  and 
Chandrasekhar  (1)  have  established  cases  for  rigorous 
kinematic  impossibility,  for  different  cases,  Herzenberg 
(2)  has  established  possibility.  Also,  axial  symmetry  of 
the  revolving  fluid  has  been  prohibited  by  the  Cowling 
theorem  (3),  with  extension  by  Backus  and  Chandra¬ 
sekhar. 

The  results  of  paleomagnetic  studies  have  strongly 
suggested  that  the  present  polarity  of  the  earth’s  main 
dipole  field  was  just  the  opposite  at  an  earlier  time  in 
geological  history.  In  fact,  polarity  reversals  are  thought 
to  have  occurred  at  various  times  during  the  existence  of 
the  earth.  This  evidence  of  polarity  reversal  from  rema¬ 
nent  rock  magnetism  is  now  generally  conceded  to  re¬ 
quire  recognition  in  any  theoretical  model  of  the  earth’s 
main  dipole  field.  Bullard  and  Gellman  (Ref.  4,  p.  258) 
and  Runcorn  (5)  have  indicated  that  the  dipole  may  be 
highly  sensitive  to  the  details  of  convection.  Thus,  polar¬ 
ity  reversal  might  conceivably  be  due  to  shifts  from  one 
steady-state  dynamo  to  another,  caused  by  extra- 
magnetofluiddynamic  changes  in  the  earth.  Despite 
much  progress,  the  steady-state  results  are  far  from 
complete. 

Though  the  theory  of  nonsteady  dynamos  is  even 
more  difficult,  many  geophysicists  believe  that  study  may 
shed  some  light  on  dipole  variation  and  polarity  reversal. 

NRL  Problem  EOI-OI;  Projects  RR  010-04-4 1-59S0  and  BuWeps  RAE 
50R-0S5/652-1/F012-1 1-001.  This  is  a  final  report  on  one  phase  of  the 
problem;  work  on  other  phases  of  the  problem  continues. 

NOTE;  One  of  the  authors  (J.H.M.)  is  formerly  of  NRL. 


A  clever  .analogy  due  to  Bullard  (6)  started  an  interest¬ 
ing  series  of  discoveries  relating  to  polarity  reversal 
capabilities  of  the  unsimplified,  nonsteady  system. 
Bullard  introduced  a  disk-dynamo  model  in  place  of  the 
vastly  more  complicated  homogeneous  dynamo  thought 
to  be  functioning  within  the  core  of  the  earth.  The 
working  principles  of  the  disk  dynamo  were  known,  and 
the  interaction  of  the  current  and  the  rotational  veloc¬ 
ity  in  this  model  could  be  precisely  formulated.  The 
disk-dynamo  model  remains  theoretical,  nevertheless, 
because  the  large  dimensions  of  the  eartl)|process  do  not 
lend  themselves  to  reproduction  on  the  small  scale 
attainable  in  the  laboratory.  The  partial  differential 
equations  of  the  magnetofluiddynamic  prototype  be¬ 
came  the  ordinary  differential  equations  of  lumped 
circuit  theory  and  rigid  rotation.  The  essential  nonline¬ 
arity  (tfxB,yx3)of  closed-system,  magnetomechani¬ 
cal  interaction  was  retained.  Admittedly,  the  disk- 
dynamo  model  has  to  be  interpf%ted  with  great  caution 
since  it  is  far  removed  from  the  details  of  the  prototype. 
With  regard  to  field  reversals,  Bullard’s  single-disk,  self¬ 
exciting,  dynamo  model  proved  somewhat  disappoint¬ 
ing,  for  he  found  that  current  oscillations  observed  in 
the  solutions  of  the  equations  were  fundamentally  incap¬ 
able  of  exhibiting  opposite  polarity.  Rikitake  (7)  worked 
out  the  theoretical  model  for  a  system  of  two,  mutually 
exciting,  disk  dynamos.  He  was  able  to  carry  out  solu¬ 
tions  of  the  equations  of  this  system  far  enough  to  estab¬ 
lish  that  the  oscillating  dynamo-currents  could  indeed 
change  sign.  Allen  (8)  published  the  solution  of  a  single 
current  for  the  same  system,  corroborating  Rikitake’s 
results.  Allan’s  solution  was  of  sufficient  length  to  reveal 
another  feature.  A  series  of  small  oscillations,  of  rip¬ 
ples,  which  remained  almost  entirely  of  one  polarity  was 
followed  by  another  series  of  opposite  polarity.  This, 
feature  seemed  to  suggest  that  the_  ripples  were  super¬ 
imposed  first  on  one  level  of  equilibrium,  and  then  on  a 
second  level  of  oppxjsite  polarity.  Thus,  in  principle,  a 
theoretical  disk-dynamo  system  model  had  been  shown 
capable  of  accounting  for  the  field  reversals  indicated  by 
the  results  of  paleomagnetic  studies.  Paleomagnetic 
evidence,  however,  seems  to  indicate  a  rather  spocial 
typje  of  variation:  intervals  of  more  or  less  random  dura- 
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tion  and  of  fairly  constant  intensity  in  one  polarity  are 
separated  by  relatively  rapid  polarity  reversals  between. 

We  have  used  an  analog  computet  to  obtain  solutions 
to  the  equations  of  the  Rikitake  two-disk  system,  with 
and  without  modification.  The  forte  of  this  computer 
is  that  continuous  solution^  of  systems  of  ordinary  non¬ 
linear  differential  equations  are  possible.  Solutions  for 
significantly  long  durations  ’  have  been  achieved  with 
relative  ease. 

The  purpose  of  this  report  is  to  show  that  a  remark¬ 
able  approximation  to  the  special  type  of  variation  is 
possible.  For  a  more  extensive  background  of  the  prob¬ 
lem  in  general,  the  reader  may  refer  to  two  recent  re¬ 
view  articles  (9,10). 

MODELS 

The  basic  Bullard-type  single-disk  dynamo  model  is 
shown  schematically  in  Fig.  1.  Voltage  is  generated 
radially  in  the  disk  as  it  rotates  in  the  coil  field  created 
by  the  current.  For  high  angular  velocities  the  smallest 
initial  field  or  current  is  sufficient  to  self-excite  the  sys¬ 
tem.  Connections  between  coil  and  disk  may  be  visual¬ 
ized  as  some  equivalent  form  of  sliding  contact.  For  the 
coil  connections  and  direction  of  rotation  shown,  this 
dynamo  is  regenerative.  If  either  is  changed  the  dynamo 
is  degenerative.  If  both  are  changed  the  dynamo  is 
again  regenerative. 


ROTATING 

DISK 


STATIONARY 

COIL 


Fig.  1  -  Bullard-type,  self-excited,  disk-dynamo  model 


Figure  2  illustrates  Rikitake’s  two-disk  system.  Disk  1 
and  Coil  I  are  coupled  electrically,  as  are  Disk  II  and 
Coil  II.  Magnetic  coupling  occurs  between  Coil  II  and 
Disk  I  and  also  between  Coil  I  and  Disk  II.  The  two 
halves,  left  and  right,  are  considered  to  be  separated 
sufficiently  to  exclude  the  presence  of  any  other  cou¬ 
pling.  Rikitake’s  dynamo  is  regenerative  for  the  coil 
connections  and  directions  of  rotation  shown. 

Figure  3(a)  shows  the  two-disk  system  of  Fig.  2  re¬ 
arranged  in  an  axial  configuration,  still  regenerative. 
For  ease  of  illustration  the  coils  have  been  drawn  closer 
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Fig.  2  -  Rikitake-type,  two-disk,  mutually  excited  dynamo 
model.  The  break  in  the  drawing  serves  to  indicate  that  the 
individual  sections,  left  and  right,  are  sufficiently  separated 
to  exclude  interaction  except  the  mutual  excitation. 
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Fig.  3  -  Axially  stacked  two-disk  system.  The  dynamo  in  (b) 
differs  from  the  dynamo  in  (a)  in  that  the  connections  of 
Coil  II  have  been  reversed. 


to  the  disks  which  are  electrically  coupled.  Figure  3(b) 
differs  trom  3(a)  m  that  connections  of  the  lower  coil 
have  been  changed  to  make  the  system  degenerative. 
The  axial  configuration  of  these  figures  requires  the 
consideration  of  additional  magnetic  coupling,  namely, 
the  mutual  inductance  between  Coil  I  and  Coil  II,  and 
the  mutual  inductances  between  Coil  I  and  Disk  I  and 
between  Coil  II  and  Disk  II. 

EQUATIONS 
Rikitake  System 

The  Rikitake  system  equations  are 

L/*  -  o/*  n* -I- R/*  =  0  (1) 

L/*  +  od*  n*  +  R/*  =  0 


(2) 
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Cft*  +  al*  I* 

-G  =0 

(3) 

ctl*  ±  aJ*  I* 

-G  =0, 

(4) 

where 


etc.,  are  time  derivatives.  I'he  asterisks  denote  that  the 
variables  are  real,  dimensional  quantities,  the  subscripts 
1  and  2  refer  to  the  two  disk-and-coil  circuits,  /  is  cur¬ 
rent,  and  fl  is  angular  velocity.  The  parameters  of  the 
system  are  symmetrical.  L  and  R  are  the  self-inductance 
and  resistance,  respectively,  of  the  circuits;  C  is  the  mo¬ 
ment  of  inertia  of  the  disks;  and  G  is  the  common,  exter¬ 
nal,  driving  torque  on  the  disks.  The  symbol  a,  called 
M  by  Rikitake  and  Bullard,  is  the  factor  which  remains 
explicit  from  the  term  iira,  which  represents  the  mutual 
inductance  between  the  coil  and  the  periphery  of  the 
disk.  For  the  two-disk  system,  this  mutual  inductance 
is  between  Coil  II  and  Disk  I  and  between  Coil  I  and 
Disk  II.  We  prefer  to  reserve  the  symbol  M  for  later 
use.  Although  somewhat  trivial,  the  lower  signs  have 
been  added  to  include  the  degenerate  type  of  disk  rota¬ 
tion  and  coil  connection. 


which  we  refer  to  as  the  A-form  and  the  P-form.  The 
A-form  relation  is: 


(9) 

“T-M 

(10) 

m 

(11) 

The  P-form  relation  is: 

i*^^Si 
'  a  Vi  ‘ 

(12) 

ft*  =  -  ft, 

<  a  ' 

(13) 

** 

ID 

(14) 

For  both  forms,  i  =  1,2. 


The  A-form  of  the  resulting  Rikitake  equations  having 
dimensionless  variables  is: 


Modified  Rikitake  System 

Generalizing  the  interaction  of  the  axial  arrangement 
of  Fig.  3  requires  the  addition  of  two  new  parameters. 
We  shall  now  let  M  represent  the  mutual  inductance 
between  the  coils.  The  factor  /3  remains  explicit  from 
the  term  27r/3,  which  represents  the  mutual  inductance 
between  Coil  I  and  the  fjeriphery  of  Disk  I,  and  between 
Coil  II  and  the  periphery  of  Disk  II.  The  equations 
become,  in  dimensional  form, 

LI*  +  MI*  -  (a/*  pi*)  n*  +  fi/f  =  0  (5) 

LI*  +  MI*  T  (oJ*  +  pi*)  n*  +  RI*  =  0  (6) 

cn*  (al*  +  pi*)  /*  -  C  =  0  (7) 

Ch*  ±  (al*  +  pi*)  I*  -G  =  0.  (8) 

The  upper  and  lower  signs  correspond  to  the  models 
in  Figs.  3(a)  and  3(b)  respectively. 

Dimensionless  Forms 

Solution  of  the  equations  by  computer  was  facilitated 
by  first  making  the  variables  dimensionless  in  order  to 
provide  for  easy  adjusttnent-of  parameters.  This  was 
done  by  means  of  the  following  Rikitake-type  relations. 


7, -7,  ft,  +  A7,  =  0 

(15) 

•  7,+7,  ft,  +  A7,  =  0 

(16) 

ft,-l-7,  7,  -1  =0 

(17) 

ft,  ±  7,  7,  -  1  =  0  , 

(18) 

where 

K-R 

The  P-form  of  the  resulting  Rikitake  equations  with 
dimensionless  variables  is: 

/, -/,n, -(-/,=  0  (19) 

/,+/,  n,-l-/,  =  0  (20) 

h, -I-/, /,  -P  =  0  (21) 

ft,  ±  /,  /,  -  P  =  0  ,  (22) 

where  P  ■  1/A*. 
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The  ^-form  of  the  modified  Rikitake  equations  for 
the  axial  system  (Fig.  3),  generalized  in  interaction,  is; 

/,  +  0/,  -  (h  +  FI,)  n^  +  KI,  =  0  (23) 

i,  +  Qi,  +  il,  +  FI,){l,  +  KI,  =  0  (24) 

h.+  (/,  +  F/.)/,  -1  =0  (25) 

h,±(/.  +  F/.)/,  -1  =0.  (26) 


(o) 


where  Q  «  MIL  and  F  ■  /3/a. 

The  /’-form  for  the  modified  system  is; 

h  +  Qi,-(i,  +  Fi,)  n,  +  /.  =  o 


},  +  Qi,zi:(I,+FI,)  a,  +  i,  =  Q 


n, +  (/,  +  /■/,)/,  -/»  =  o 

h,±(/.  +  F/,)/,  -p  =  o. 

RESULTS 


(30) 


Fig.  4  -  (a)  Two  individual  currents  I,  and  /,  for  P  =  1.000. 
Initial  conditions;  I,  =  3.000,  /,  =  1.914,  fl,  =  4.000,  fl,  =  0.037. 
(b)  Repeat  runs  of  the  same  curves,  and  their  sum  +  /,. 


The  study  performed  on  the  analog  computer  for  a 
variety  of  parameters  and  initial  conditions  resulted  in 
many  solutions  of  paleomagnetic  interest.  In  the  search 
for  interesting  behavior,  solutions  were  observed  after 
each  change  of  a  fairly  systematic  variation  of  param¬ 
eters  and  initial  conditions.  Occasionally,  rapid  changes 
of  style  were  observed  within  narrow  limits  of  param¬ 
eter  and  initial  condition  changes.  Although  good 
coverage  is  thought  to  have  been  obtained,  the  search 
was  not  automated  and  cannot  be  claimed  to  have  been 
completely  exhaustive.  The  recorded  results  were  ob¬ 
tained  in  terms  of  current  versus  time  (both  dimen¬ 
sionless)  where  current  represents  the  field. 

Solutions  to  the  Rikitake  Equations  (P-Form) 

The  first  group  of  solutions  are  from  the  dimension¬ 
less  P-form  equations  of  the  Rikitake  system,  namely, 
Eqs.  (19)-(22),  using  the  upper  signs  (regenerative  case) 
and  the  value  P  =  1 .000.  Figure  4(a)  shows  the  individual 
currents  in  a  typical  reversing  type  of  behavior.  As  can 
be  seen,  both  curves  consist  of  series  of  small  oscilla¬ 
tions  in  one  polarity  followed  by  another  series  usually 
in  the  opposite  polarity.  Thus  piolarity  reversals  have 
occurred  between  the  series  of  small  oscillations.  The 
upper  curve  (/,)  resembles  strongly  the  one  given  by 
Allan  (8).  However,  for  this  curve,  portions  of  the 


small  oscillations  in  a  series  of  one  polarity  frequently 
extend  across  the  zero<urrent  line  into  the  opposite 
polarity.  Sig^nificantly,  for  the  current  in  the  lower 
curve  (/,),  the  series  of  small  oscillations  change  polar¬ 
ity  in  their  entirety.  The  small  oscillations,  or  ripples, 
are  distinctly  separated  from  the  zero  line  until  reversal 
action  takes  place,  i.e.,  the  small  oscillations  do  not  cross 
the  zero  line.  We  consider  this  distinctly  separated  typ>e 
of  solution  to  be  for  more  appropriate  for  representing 
the  rather  special  typo  of  earth’s  field  reversal  behavior 
suggested  by  paleomagnetic  data.  Figure  4(b)  shows 
three  solutions.  Two  of  these  are  repoat  runs  of  the 
individual  currents  for  comparison  with  the  third,  which 
is  the  sum  of  the  two.  Several  attempts  were  found 
necessary  to  sequentially  record  this  series  to  obtain 
results  which  would  be  equivalent  to  a  simultaneous 
recording  of  the  three  curves.  Even  so,  a  slight  scatter 
is  apparent  beyond  the  value  of  60  on  the  dimension¬ 
less  time  scale. 

The  highly  sensitive  nature  of  the  equations  in 
regions  of  zero  current  is  illustrated  in  Fig.  5(a)  for 
the  parameter  value  P  =  1.  The  second  run  departs 
drastically  from  the  first  at  t  =  48.  The  departure  is 
seen  to  consist  of  opposite-polarity  reversals  which 
approximate  a  mirror  image  of  the  first  run.  The  small 
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Fig.  5  -  Two  runs  of  the  current  sum  /,  +  /,  for  (a)  P  —  1.000 
and  (b)  P  =  0.550.  Initial  conditions:  /,  =  3.000,  /,  =  1.914, 
n,  =  4.000,  and  n,  =  0.037. 

difference  which  is  observable  at  t  =  35  may  have  been 
the  cause  of  the  drastic  departure  at  t  =  48.  The  highly 
sensitive  regions  at  t  =  9  and  22  were  successfully 
repeated  on  the  second  computer  run.  The  single 
excursion  above  the  zero  line  at  t  =  9  is  considered  to 
be  an  initial  reversal  to  positive  polarity  followed  by 
an  immediate  return  to  negative  polarity. 

Although  repeatability  varied  somewhat  with  param¬ 
eters  and  initial  conditions,  these  curves  are  rather 
typical.  That  is,  approximately  three  or  four  high- 
sensitivity  regions  is  the  limit  that  can  usually  be  suc¬ 
cessfully  retraced.  Certainly  the  solutions  are  subject 
to  reversal  whenever  they  enter  the  highly  sensitive 
region  about  zero  current.  Somewhat  longer  repeata¬ 
bility  is  illustrated  in  Fig.  5(b)  for  a  lower  value  of  the 
parameter  P.  Here,  however,  just  one  highly  sensitive 
region  has  been  encountered,  exclusive  of  the  polarity 
change  at  the  very  start. 

Although  several  individual<urrent  solution  record¬ 
ings  are  presented,  the  sum  of  the  currents  was  recorded 
for  many  of  the  solutions  in  this  first  group.  One  rea¬ 
son  was  that  in  the  early  phases  of  the  work  we  attached 
significance  to  the»  sum  for  axially  stacked  disk  and 
coil  models.  Later  we  realized  that  the  sum  in  an  axial 
variation  of  the  Rikitake  system  would  not  be  very  real¬ 


istic  without  inclusion  of  the  resulting  interaction. 
There  was,  however,  a  practical  advantage.  The  sum 
of  the  currents  in  the  Rikitake  system,  for  any  con¬ 
ditions  observed,  always  displayed  characteristics  sim¬ 
ilar  to  those  of  the  individual  current  which  is  distinctly 
separated  from  zero.  This  is  due  to  the  small  contribu¬ 
tion  from  the  other  individual  current.  A  distinctly 
separated  curve  appeared  on  /,  for  some  runs  and  on 
/j  for  others,  depending  on  the  particular  parameters 
and  initial  conditions  chosen.  Thus  the  sum  provided 
a  quick  monitor  of  the  most  interesting,  distinctly  sep¬ 
arated,  individual  current  when  utilizing  a  single 
recorder. 

The  effects  of  P-parameter  variation  for  one  set  ot 
initial  conditions  may  be  seen  from  comparison  of 
Figs.  6(a),  (b),  and  (c).  The  relatively  low  P-value  of 
Fig.  6(a)  results  in  there  being  many  more  ripples  dur¬ 
ing  the  intervals  between  reversals.  Also,  the  current 
amplitude  is  smaller.  The  nature  of  the  ripples  of 
Fig.  6(a)  may  be  seen  more  clearly  by  referring  to  Fig. 
5(b),  both  for  a  P-value  of  0.550.  Although  the  solu¬ 
tion  of  Fig.  7  was  based  up>on  a  slightly  different  initial 
condition  (ft,  =  0.0),  it  serves  to  show  typical  detail  of 
the  ripple  content  of  Fig.  6(b).  Both  solutions  are  for 
a  P-value  of  1.000.  P-values  lower  than  about  0.550 
resulted  in  solutions  (not  shown)  which  no  longer 
reversed. 

Figure  8(a)  illustrates  an  interesting  type  of  behav¬ 
ior.  The  individual  current  shows  preference  for 
extended  intervals  in  one  polarity  separated  by  brief 
multiple,  or  even  single,  excursions  into  the  opposite 
polarity.  With  some  changes  in  the  initial  conditions, 
the  individual  current  of  Fig.  8(b)  exhibits,  initially, 
similar  behavior  to  Fig.  8(a),  but  in  the  opposite  polar¬ 
ity.  It  then  reverses  polarity  and  settles  to  a  very  long 
oscillatory  mode.  We  have  sometimes  visually  monitored 
solutions  of  this  tyjje,  long  after  the  recorder  came  to 
the  end  of  its  sheet  of  paper,  by  using  digital  and  elec¬ 
tronic  voltmeters.  Nothing  was  observed  which  might 
indicate  that  the  oscillatory  mode  would  ever  be  dis¬ 
turbed.  Eventually,  dnft  and  other  signs  of  deterio¬ 
ration  of  computer  performance  became  evident,  as 
could  be  expected  from  the  outrageously  long  comput¬ 
ing  time  demand. 

Figure  9(a)  shows  a  solution  of  special  paleomagnetic 
interest.  Although  the  time  scale  may  not  be  that 
thought  precisely  desirable  from  rock  magnetism 
studies,  the  run  shows  that  quite  ideal  representations 
are  indeed  possible  with  a  disk-dynamo  system  model. 
Figures  9(b)  and  (c),  obtained  on  a  later  date  for  the 
same  conditions  used  in  Fig.  9(a),  illustrate  another 
point  of  the  repeatability  difficulties.  When  repeating 
the  longer  runs,  the  sequence  of  reversals  changed 
in  somewhat  random  fashion,  Nevertheless,  the  over¬ 
all  behavior  patterns  associated  with  the  conditions 
used  can  usually  be  observed  in  the  repeated  solution. 
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Fig.  6  -  Current  sum  /,  +  /,  for  (a)  P  =  0.550,  (b)  P  “  1.000, 
and  (c)  P  =  1.400.  Initial  conditions:  /,  ”  3.000,  /,  =  1.914, 
n,  =  4.000,  and  O,  =  0.037. 


Fig.  7  -  Expanded  scale  of  current  sum  for  P  =  1.000  (see 
Fig.  6(b))  to  illustrate  typical  ripple  content.  Initial  condi¬ 
tions:  /,  =  3.000,  /,  =  1.914,  ft,  =  4.000,  and  ft,  =  0.000  (sUghtly 
different  value  than  that  of  Fig.  6(b)). 
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Fig.  8  -  Individual  current  for  P  =  0.580.  (a)  Initial  conditions: 
/.  =  3.000,  /,  =  1.914,  ft,  =  4.000,  and  ft,  =  0.037;  (b)  Initial 
conditions:  /,  =  3.000,  7,  =  1.915,  ft,  •  3.802,  and  ft,  =  0.0. 


The  pattern  of  the  more  frequent  reversals  toward 
the  right-hand  side  of  these  curves  seemed  quite  de¬ 
pendent  upon  the  exact  manner  in  which  the  first 
reversal  to  positive  polarity  occurred.  Probably  we 
were  somewhat  fortunate  to  capture  the  long  stay 
above  the  line  in  Fig.  9(a).  On  further  attempts  (not 
shown)  we  were  able  to  reproduce  Fig.  9(a)  as  far  as 
t  =  675.  However,  small  refinements  of  the  earlier 
settings  proved  necessary  to  accomplish  this. 
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Fig.  9  -  Three  separate  runs  for  0.715  with  initial  conditions: 
/,  =  3.300,  /,  =  0.0,  n,  =  0.0,  and  O,  =  2.400.  (a)  illustrates  a 
long  duration  in  the  upper  polarity  which  occurred  after  three 
*  rapid-order  reversals;  (b)  and  (c)  illustrate  the  general  simi¬ 
larity  of  style  for  the  more  random  reversab  of  the  longer 
runs,  (a)  and  (c)  are  for  the  current  sum  /,  +  /,,  and  (b)  is 
for  the  individual  current  /,. 


Figures  10(a)  and  10(b)  present  the  only  solutions 
appearing  in  this  report  which  were  recorded  simul¬ 
taneously  during  a  single  run,  using  two  recorders. 
A  comparison  of  these  figures  illustrates  our  obser¬ 
vations  concerning  the  relation  between  the  individual 
currents  and  their  sums  for  the  Rikitake  system,  i.e., 
the  sum  of  the  currents  was  found  representative  of 
one  of  the  individual  currents  for  all  of  the  param¬ 
eters  and  initial  conditions  attempted.  The  shorter 
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Fig.  10  -  Simultaneous  recording  of  (a)  the  individual  current 
/,  and  (b)  the  current  sum  /,  -F  /,,  both  for  P  —  0.600  and  the 
initial  conditions:  /,  =  3.000,  /,  ==  1.914,  O,  =  4.000,  and  f), 
=  0.037. 

• 

length  of  run  recorded  for  the  current  sum  in  Fig. 
10(b)  is  a  result  of  the  smaller  size  of  the  second 
recorder. 

Solutions  to  the  Modified  Equations  (A(-Fonn) 

The  next  group  of  curves  presented  in  this  reptort 
are  recorded  solutions  of  Eqs.  (23)-(26)  using  the  upper 
signs.  These  are  the  dimensionless  ^(-form  of  the  equa¬ 
tions  which  contain  the  generalized  interaction  terms 
suggested  by  the  axial  system  shown  in  Fig.  3(a).  Only 
the  sum  of  the  currents  was  recorded  for  this  second 
group  of  solutions.  Here,  however,  the  sum  is  not 
necessarily  an  indicator  of  the  behavior  of  either  of 
the  individual  currents.  With  the  addition  of  mutual 
coupling,  the  current  sum  or  the  resultant  field  of  all 
the  currents  in  the  system  seems  to  be  the  only  con¬ 
sistent  representation  of  the  total  field.  Figures  11(a) 
and  (b)  illustrate  just  one  of  the  many  interesting  types 
of  behavior  encountered  in  the  second  group  of  solu¬ 
tions.  The  A-value  of  zero  for  Fig.  11(a)  is  probably 
unrealistic.  A  changing  frequency  of  the  small  oscil¬ 
lations  may  be  observed  in  both  of  these  solutions. 
Figures  12(a)-(f)  illustrate  a  variety  of  styles  of  behav¬ 
ior  found  in  these  solutions.  It  is  apparent  that  the 
particular  set  of  initial  conditions  and  parameters  does 
not  forever  determine  the  reversal  pattern. 
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Fig.  1 1  -  The  current  .sum  /,  +  /,  for  (a)  K  =  0.0,  F  =  0.0075, 
and  Q  =  0.0350,  and  (b)  K  =  0.0025,  F  =  0.o875,  and  0  =  0.0850. 
Initial  conditions:  /,  =  3.000,  /,  =  1.914,  fl,  =  4.000,  and  fl,  =  0.0. 


Much  larger  values  of  the  F-  and  ^-parameters  re¬ 
sulted  in  solutions  (not  shown)  which  did  not  reverse, 
except  for  an  cKcasional  polarity  change  at  the  very 
start.  The  ripples  usually  settled  to  constant,  non¬ 
zero,  steady-state  values. 

Some  solutions  (not  shown)  were  obtained  for  the 
model  of  Fig.  3(b)  using  Eqs.  (27)-(30)  with  the  lower 
signs  (degenerative  case).  In  general,  these  solutions 
settled  rapidly  to  zero  from  the  initial  conditions.  One 
band  of  initial  conditions  and  parameters  was  found 
which  resulted,  at  first,  in  a  reversal  with  almost  every 
oscillation,  but  then  soon  settled  to  zero. 

DISCUSSION 

Attempts  to  predict  in  advance  the  time  of  occur¬ 
rence  of  the  next  reversal  in  the  solutions  proved  un¬ 
successful.  While  increasing  ripple  amplitude  of  the 
currents  is  frequently  evident  prior  to  reversal,  no 
reliable  criteria  could  be  established.  A  linearization 
technique  was  tried  with  the  hope  that  differences 
between  linear  and  nonlinear  solutions  might  fitretell 
the  (Kcurrence  of  the  next  reversal.  However,  no  ad¬ 
vanced  indication  was  obtained  from  this  technique. 


Figure  13  is  an  idealized  pictorial  representation 
of  some  of  the  important  characteristics  of  the  earth’s 
field  reversals  as  indicated  from  paleomagnetic  studies. 
Originally,  it  seemed  conceivable  in  the  present  study 
that  parameter  and  initial  condition  values  might  be 
found  which  would  give  solutions  roughly  constant  in 
one  polarity,  between  the  reversals.  Th^  solid  line  of 
the  sketch  of  Fig.  13  illustrates  this  idea.  However,  no 
recurring-reversal  solutions  without  ripple  were  found, 
even  though  an  extensive  parameter  and  initial  con¬ 
dition  search  was  made.  Actually,  the  presence  of  rip¬ 
ples  during  the  interval  in  one  polarity,  illustrated  by 
the  dashed  line  of  the  sketch,  does  not  clash  with  paleo- 
magnetic  evidence.  Allan  (8)  cites  several  references 
containing  evidence  of  field  intensity  variation.  A  lack 
of  additional  evidence  may  be  due  to  inability  of  paleo¬ 
magnetic  techniques  to  resolve  the  ripple  content. 

The  interval  of  time  during  reversal,  from  paleo¬ 
magnetic  evidence,  should  be  of  the  order  of  10,000 
years.  The  field  should  remain  in  one  polarity  for  in¬ 
tervals  of  100,000  years  to  2  or  3  million  years,  with 
500,000  years  considered  typical  (10,11).  For  compar¬ 
ison,  the  time  scales  of  the  solutions  presented  may 
be  changed  to  dimensional  units  (seconds)  using  Eq. 
(14)  for  the  P-form  and  Eq.  (11)  for  the  K-form.  It  is 
difficult  to  estimate  reliable  values  of  the  quantities 
needed  for  these  time-scale  conversions,  particularly 
of  L  and  G.  If  Eq.  (11)  is  multiplied  by  l/K,  there 
results  the  optional  form 


'•-C©..  (51) 

Thus  the  time-scale  conversion  for  the  K-iorm  solu¬ 
tions  differs  from  the  P-form  conversion  only  by  the 
factor  K.  The  alternate  comparison  may  be  obtained 
by  multiplying  Eq.  (14)  by  1/VP,  giving  the  optional 
form 


t*  =  VP^t.  (32) 

It  is  thought  that  the  task  of  estimating  reliable  values 
of  L  and  R,  or  of  C,  L,  G,  and  o,  for  representing  the 
earth’s  core  is  best  left  to  geophysicists  specializing  in 
these  matters.  Hence,  no  direct  attempt  will  be  made  to 
convert  the  solutions  to  dimensional  time.  Nevertheless, 
it  is  worth  pointing  out  that  the  ripple  frequencies  of 
the  first  group  of  solutions  (Figs.  4-11)  are  within  a 
factor  of  2  or  3  of  those  of  Rikitake’s  and  Allan’s.  Thus 
it  may  be  inferred  that  the  dimensional  times  of  these 
solutions  to  the  Rikitake  equations  would  be  con¬ 
sidered  insufficient  by  a  factor  of  about  500  for  repre¬ 
senting  paleomagnetic  times. 

The  second  group  of  solutions  correctly  represents 
the  behavior  of  a  two-disk  dynamo  system  with  mutual 
interaction  in  axial  arrangement,  and  they  may  be  of 
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Fig.  12  -  The  current  sum  /,  +  /,  for  (a)  K  =  0.0002,  F  =  0.0400,  and  0  =  0.4000,  with  initial  condi¬ 
tions:  4  =  6.6000,4  =  0.0040,  01  =  0.8000,  and  =  0.8000;  (b)  A  =  1.5  x  10-‘,f =0.0210,  and  0  = 
0.4000,  with  initial  conditions:  /■  =  6.1800,  It  =  2.3800,  Oi  =  O.SOOO,  and  Oi=  1.9400;  (c)X=0.0001, 
F  =  0.0400,  and  Q  =  0.4000,  with  initial  conditions:  A  =5.6000, /i  =  4.4000, 01=3.4400,  and  ns  =  0.0; 
(d)  K  =  0.0001,  F  —  0.0200,  and  Q  =  0.4000,  with  initial  conditions:  Ii=  1. 4000, A” 0.7660, Oi= 
1.8000,  and  Oi  =  0.0;  (e)  K  =1X  lO”’,  F  =  0.0285,  and  0=0.5000,  with  initial  conditions:  A  =  0.7640, 
4  =  1.2000,  n,  =  0.0040,  and  n,=  1.6000;  andj(f)X- 1  X  10-»,F= 0.1050,  andO-0.4000,  with  initial 
conditions;  A  ~  2.7880,  A  =  4.2380,  tli  =  1.8080,  and  tU  =  1.0440. 
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Kig.  13  -  Idealized  pictorial  representation  of  important  characteristics  of  earth’s  field 
reversals  indicated  from  paleomagnetic  studies 


some  interest  ptr  se.  However,  Backus  (12)  has  found 
that  even  nonsteady  axisymmetric  dynamos  are  incap¬ 
able  of  long-term  support  of  external  dipole  moments. 
Obviously,  the  axial  arrangement  of  the  dynamo  sys¬ 
tems  of  Fig.  3  must  not  be  taken  too  literally.  The 
interactions  must  be  construed  in  a  highly  schematic 
way,  not  confined  to  the  configuration  which  suggested 
them. 

Paleomagnetically  interesting  behavior  of  the  mod¬ 
ified  equations  was  obtained  from  the  /C-form  using 
relatively  low  values  of  K.  Those  solutions  with  very 
low  values  of  K  are  probably  the  most  realistic  for 
representing  the  earth’s  core  (see  Ref.  (7),  p.  104). 
However,  the  dimensionless  ripple  frequencies  of  the 
second  group  of  solutions  are  not  greatly  different 
from  those  obtained  for  the  unmodified  equations. 
If  Eq.  (31)  were  used,  these  very  low  K-values  would 
cause  the  dimensional  times  of  the  modified  solutions 
to  be  unreasonably  small  for  interpretation  of  paleo¬ 
magnetic  behavior. 

Probably  further  efforts  should  be  made  to  reevaluate 
the  parameters  L,  R,  a,  C,  and  G.  Should  significant 
revision  result,  a  more  careful  search  for  interesting 
solution  behavior  for  the  new  values  would  seem  to  be 
in  order. 

Despite  the  difficulties  inherent  in  determining  proper 
dimensional  time  scales,  many  of  the  solutions  presented 
show  remarkable  similarity  to  the  reversals  of  the 
earth's  main  dipole  field.  Clearly  there  is  a  choice  of 


an  extremely  wide  range  of  intervals  in  one  polarity 
regardless  of  any  final  adjustment  to  the  time  scales. 
Acceptable  ratios  (10/1  to  100/1)  of  the  interval  in  one 
polarity  to  the  interval  during  reversal  are  found. 
Paleomagnetic  evidence  further  suggests  that  single 
reversals  or  excursions  of  short  duration  into  opposite 
polarity  may  have  occurred.  This  feature  is  present  in 
most  of  these  solutions. 

The  use  of  an  analog  computer  has  resulted  in  a 
significant  extension  of  disk-dynamo  model  solutions. 
In  repeat  runs  using  the  same  parameter  values, 
exact  reproducibility  was  not  obtained  beyond  several 
ripples  of  the  solutions.  Presumably  this  is  due  to 
rep>eated  approaches  of  the  solutions  to  a  mathematical 
singularity.  Increased  accuracy  and  reproducibility, 
resulting  in  a  postponement  of  scatter,  could  be  ex¬ 
pected  from  a  digital  computer  utilizing  a  large  number 
of  digits.  But  the  analog  method  insures  that  the  exact 
mathematical  formulation  does  not  mislead  one  into 
exp>ecting  unnatural  physical  exactness. 

The  simplification  from  the  homogeneous  dynamo, 
thought  to  exist  in  the  earth’s  core,  to  disk-dynamo 
models  is  a  drastic  one,  so  an  interpretation  of  the 
results  of  the  simplified  models  must  be  cautious.  How¬ 
ever,  the  solutions  presented  in  this  report  emphasize 
the  capability  of  disk-dynamo  system  models  to  per¬ 
form  reversals  satisfying  paleomagnetic  criteria.  Thus 
the  magnetofluiddynamic  origin  of  reversals  of  the 
earth’s  field  becomes  more  plausible. 
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